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We discuss the equation of state for 2 flavor QCD at non-zero temperature and den- 
sity. Derivatives of In Z with respect to quark chemical potential fi q up to fourth 
order are calculated, enabling estimates of the pressure, quark number density and 
associated susceptibilities as functions of n q via a Taylor series expansion. It is 
found that the fluctuations in the quark number density increase in the vicinity of 
the phase transition temperature and the susceptibilities start to develop a pro- 
nounced peak as fi q is increased. This suggests the presence of a critical endpoint 
in the (T, fi q ) plane. 
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1. Introduction 
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Remarkable progress for QCD thermodynamics has been recently made by 
numerical studies of lattice QCD with small but non-zero baryon density 1 . 
The pseudocritical line has been investigated in the low density regime 2 ' 3 ' 4 . 
The phase transition for 2 flavor QCD is known to be crossover at /j, q = 
and expected to become a first order phase transition at a critical endpoint, 
and it might be possible to detect the endpoint experimentally via event- 
by-event fluctuations in heavy-ion collision experiments. 

In this report, we discuss the equation of state at non-zero baryon num- 
ber density. The study of the equation of state gives the most basic infor- 
mation for the experiments. Quantitative calculations of thermodynamic 
quantities such as pressure and energy density are indispensable. In par- 
ticular, since the number density fluctuation should be large around the 
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endpoint 5 , the susceptibility of the quark number density, given by the 
second derivative of pressure with respect to ji q , is an important quantity 6 . 

Several studies of the quark number susceptibility have been performed 
at \i q — 7 . Moreover, measurements of pressure and energy density at 
H q ^ were done by Ref. 8 using the reweighting method, which allow 
the investigation of thermodynamic properties at non-zero baryon density 
This approach, however, does not work for large \i q and large lattice size 
due to the sign problem. 

Here, we adopt the following strategy. We compute the derivatives of 
physical quantities with respect to \i q at [i q = 0, and determine the Taylor 
expansion coefficients in /i 9 3 ' 9 , in which the sign problem does not arise. 
Because the pressure is an even function of /j, q , the /x^-term is leading and 
/i q -term is the next to leading, and in fact only these two terms are non-zero 
in the high temperature limit. We compute the Taylor expansion coeffi- 
cients up to fourth order. Using the second derivative of pressure (energy 
density), we investigate the relation between the line of constant pressure 
(energy density) and the phase transition line 3 . The fourth order term en- 
ables us to evaluate the /Xq-dependence of the quark number susceptibility 
near fi q = . By estimating the change of the susceptibility, we discuss the 
possibility of the existence of the critical end point in the phase diagram of 
T and n q w . 

2. Taylor expansion in fi q 

Pressure is given in terms of the grand partition function Z(T, V, /i q ) by 
p/T 4 = (l/UT 3 )lnZ. However, the direct calculation of InZ is difficult, 
hence most of the work done at [i q = for the calculation of pressure is done 
by using the integral method 11,12 , where the first derivative of pressure is 
computed by simulations, and the pressure is obtained by integration along 
a suitable integral path. For \i q ^ 0, direct Monte Carlo simulation is not 
applicable; in this case we proceed by computing higher order derivatives 
of pressure with respect to fi q /T at fi q = 0, and then estimate p(fJ, q ) using 
a Taylor expansion, 



where c„ = (l/n!)d n {p/T 4 )/d(fj, q /T) n \ flq=0 . These derivatives can be com- 
puted by the random noise method, which saves CPU time, and also it can 
be proved that the odd terms are exactly zero. Furthermore, we do not 
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Figure 1. Coefficients of Taylor expansion, C2 (left) and C4 (right). Tq is T c at fj, q = 0. 



need simulations at (T, fjb q ) = (0, 0) for the subtraction to normalize the 
value of p, since the derivatives of p| T=0i([t(;=u with respect to \x q are, of 
course, zero. This also reduces CPU time. 

The quark number density n q and quark number susceptibility Xq 
are obtained by the derivatives of pressure; n q /T 3 = d(p/T 4 )/d(^i q /T), 
X q /T 2 = d 2 (p/T 4 )/d(p q /T) 2 . We compute the pressure up to 0{n\) using 
2 flavors of p4-improved staggered fermion 13 at ma = 0.1 on a 16 3 x 4 
lattice. Then the quark number density and quark number susceptibility 
are obtained up to 0(p? q ) and 0(/j 2 ), respectively. The details of the simu- 
lations are given in Ref. 10. In Fig. 1, we plot the data for C2 (left) and c 4 
(right). Both of them are very small at low temperature and approach the 
Stefan-Boltzmann (SB) limit in the high temperature limit, as we expected. 
The remarkable point is a strong peak of C4 around T c . 

We can understand this peak via two arguments. One is a prediction 
from the hadron resonance gas model 14 , which is an effective model of the 
free hadron gas in the low temperature phase. The model study predicts 
C4/C2 = 0.75 and our results are consistent with this prediction for T <T C ; 
in fact, as T increases C4/C2 remains constant until T w T c , whereupon it 
approaches the SB limit. 

The other point is from a discussion of the convergence radius of the Tay- 
lor expansion. We expect that the crossover transition at [i q — changes 
to first order transition at a point n q /T c ~ O(l) 2 ' 15 . Then, the analysis 
by Taylor expansion must break down in that regime, i.e. the convergence 
radius should be smaller than the value of \x q jT at the critical endpoint. 
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We define estimates for the convergence radius by p n = y/\c n /c n+ 2\- We 
compute po and pi from cq = p/T 4 (Q), c 2 and C4. It is found that both po 
and P2 are quite large at high temperature as expected from the SB limit, 
pf B ~ 2.01, pf B ~ 4.44. On the other hand, around T c , these are O(l), 
since c 2 and c 4 are of the same order, so that our results around T c suggest 
a singular point in the neighborhood of p q /T c = 1. 

3. Equation of state at fi q ^ 

Next, we calculate pressure and quark number susceptibility in a range 
of < fJ. q /T < 1 which is within the radius of convergence discussed 
above, using the data of c 2 and c 4 ; A(p/T 4 ) = p(T, p q )/T 4 - p(T, 0)/T 4 = 
c 2 (p q /T) 2 + Ci {p q /T) 4 + 0(jj«), and Xq /T 2 = 2c 2 + l2 Ci (p q /T) 2 + 0(p 4 ). 

We draw A(p/T 4 ) for each fixed p q /T in Fig. 2 (left) and find that 
the difference from p\^ q= o is very small in the interesting regime for heavy- 
ion collisions, p q /T sa 0.1 (RHIC) and p q /T « 0.5 (SPS), in comparison 
with the value at \i q = 0, e.g. the SB value for 2 flavor QCD at p q = 
0: p SB /T 4 ~ 4.06. The effect of non-zero quark density on pressure at 
p q /T = 0.1 is only 1%. Also, the result is qualitatively consistent with that 
of Ref. 8 obtained by the reweighting method. 

Moreover, together with the data of derivative with respect to /3, we 
can discuss the lines of constant pressure [energy density] at T c . From 
dp(T, p 2 ) = 0, the slope at p q = is 

dT d(p/T 4 ) I ( d{p/T 4 ) Ap\ 

<W " d{p? q ) / I W + W' 

We obtain T 2 d 2 (p/T 4 )/dp 2 q = 0.693(5), and use the data in Ref. 11, p/T 4 = 
0.27(5), Td(p/T 4 )/0T = -[a-^da/d/S^idip/T^/dlS) = 2.2(6), at C 
for ma = 0.1. The same calculation is also performed for the energy density. 
We find that the slope of the constant pressure [energy density] line is 
T{AT/A{p 2 q )) = -0.107(22)[-0.087(23)]. Since the slope of T c in terms of 
p\ is T c (dT c /d(p 2 q )) = -0.07(3) 3 , this result suggests that the line of 
constant pressure or energy density is parallel with the phase transition 
line, and the ^-dependence of the pressure or energy density at T c (p q ) is 
very small 3 . 

Figure 2 (right) is the result for x q /T 2 at fixed p q /T. We find a pro- 
nounced peak for m q /T > 0.5, whereas \q does not have a peak for p q = 0. 
This suggests the presence of a critical endpoint in the (T,p q ) plane. 

This discussion can be easily extended to the charge fluctuation \c- 
The spike of \c at T c is weaker than that of x q , which means the increase 
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Figure 2. Difference of pressure from /i q = (left) and Quark number susceptibility 
(right) as a function of T for each fixed /J, q /T. To is T c at fj, q = 0. 

of the charge fluctuation is smaller than that of the number fluctuation as 
H q increases 10 . 
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